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Abstract 

An ultracold gas is considered, loaded into a lattice, each site of which is formed by 
a double-well potential. Initial conditions, after the loading, correspond to a nonequi- 
librium state. The nonlinear dynamics of the system, starting with a nonequilibrium 
state, is analysed in the local-field approximation. The importance of taking into ac- 
count attenuation, caused by particle collisions, is emphasized. The presence of this 
attenuation dramatically influences the system dynamics. 
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1 Introduction 



Present-day experiments provide the possibility of creating a large variety of optical and 
magnetic lattices with loaded ultracold Bose and Fermi gases [1-6]. One of the main advan- 
tages of such systems is the potentiahties of their almost complete control. The parameters 
of the lattices themselves can be varied in a wide range, as well as there exists a great di- 
versity of characteristics of ultracold gases, both Bose [7-13] as well as Fermi [14-16]. Such 
a high level of control suggests the feasibility of numerous applications. For many of these 
applications, it is important to have a correct description of the system dynamic properties. 

Recently, the variety of lattices has been expanded with the experimental realization of 
double-well optical lattices, each site of such lattices being formed by a double-well optical 
potential [17]. Such more complicated double- well lattices possess a rich variety of ther- 
modynamic equilibrium states [18]. Dynamics of atoms in these complex lattices should 
also be more rich, presenting novel possibilities that could be employed for various practical 
purposes, for instance, in quantum information processing and quantum computing. 

Atomic dynamics inside a separate double- well have been studied earlier [19-21]. The 
aim of the present paper is to analyze dynamics not inside a separate double- well, but in the 
whole double-well lattice, composed of many double-well sites. 

We consider atoms in a double-well lattice, which in equilibrium would form the Mott 
insulating state. But initially atoms are loaded in a nonequilibrium state. Describing their 
dynamics, it is important to invoke realistic approximations, taking into account attenuation 
caused by atomic collisions. Resorting to a simple mean-field approximation, usually em- 
ployed for treating atomic dynamics inside a separate double well, would result in perpetual 
oscillations, never coming to a stationary state corresponding to equilibrium. Therefore, 
including attenuation is of crucial importance. 

The plan of the present paper is as follows. First, wc derive the effective Hamiltonian 
for ultracold gases in an insulating double- well lattice (Section 2). Then we consider the 
nonlinear evolution equations, taking account of particle interactions and attenuation, for 
the dynamic variables describing particle tunneling, Josephson current, and atomic displace- 
ments (Section 3). In Section 4, we investigate the existence and stability of attractors, and 
present the results of numerical calculations. Finally, Section 5 is conclusion. 



2 Effective Hamiltonian 

The starting point is the standard energy Hamiltonian in the Heisenberg field representation, 

H = J ^t(r)i7^(r)^(r) dr + ^ J ^^(r)^^(r')$(r - r')V'(r')^(r) drdr' , (1) 

in which il}{r) is a field operator, where the time dependence, for brevity, is omitted, but 
keeping in mind that ip{r) = '0(r, is the Heisenberg field operator; $(r) = ^(~r) is an 
interaction potential; and 

where the lattice potential Vl(i') is assumed to have the double- well structure in each lattice 
site marked by a lattice vector cij. For generality, we do not specify a particular form of the 
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interaction potential $(r), since in experiments it can have different forms. It can be a local 
potential [7-11], or a long-range dipolar potential, or a combination of a contact and dipolar 
potentials [22]. The strength of atomic interactions can also be varied in a very wide range 
by using Feschbach resonance techniques [23,24]. 

The lattice is assumed to be sufficiently deep in order to realize the Mott insulating state, 
when the wave functions </?„(r — eij), which are the solutions of the eigenproblem 

HL{r - ai)(pn{r - a^) = E^ip^ir - a^) , (3) 

are well localized around the related lattice sites a^. Then the jumps of atoms between 
different sites are suppressed, as it should be for an insulating state. However the tunneling 
between the wells of a double-well potential, associated to each lattice site, can be arbitrary 
and can be regulated by varying the shape of the double well [25]. 
Expanding the field operator as 

i'i^) = Cm<^n(r - 3ii) , (4) 
in 

we can transform Hamiltonian (1) to the representation of operators q^, whose commutation 
relations are prescribed by the sort of atoms, which can be either bosons or fermions. In 
an insulating state, the filling factor is to be integer, which we set to one. This imphes the 
homeopolarity conditions 

^ ^ C-in^in — 1 ) QmCin — . (5) 

n 

Substituting expansion (4) into Hamiltonian (1) and retaining overlaps up to second order, 
gives 

rV _ \ ' p „ I ^ \ ' \ ' ^rmnm'n' J J /p\ 
-'^n^in^in ' „ 2-^ ^ ij ^im'^jn'^jm"^in' i ) 



where 
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i^j mnm'n' 



irmnm n — ^mnm n i ^mnn 
ij ~ ijji ijij 



is expressed through the matrix elements of the interaction potential with respect to the 
wave functions V9n(i" — cij). 

The system temperature is assumed to be close to zero, so that, of the double-well 
spectrum, only two lowest energy levels may be taken into account. These levels, Ei and 
E2, are enumerated with the index n — 1 and n — 2. The wave functions of the lowest 
levels, corresponding to a double-well potential, are such that the ground-state function is 
symmetric with respect to spatial inversion, while the first excited state is antisymmetric 
[21], 

V^i(-r) = <^i(r) , <^2(-r) = -<^2(r) . (7) 

Both functions (pi{r) and (r) are real. 
By introducing the pseudospin operators 

St^l (4iQ2 + 4ca) , (8) 
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we have 



cliCii — + , c|2Ci2 — 2 ~ ^ 



4,c^ = s^-tsy, 4,c,^ = + 1 sf . (9) 

In order to better understand the meaning of the pseudospin operators (8), we may define 
the left, CiL, and the right, cm, operators by the relations 

Cil = {CiL + Ciu) , Ci2 = -y^ {CiL - Ciu) . (10) 

Then the operators (8) acquire the form 

{ciLdL - clRCiR) ■ (11) 

^Prom here, we see that Sf describes the tunnehng intensity between the left and right wells 
of a double-well potential located at the i-site. The operator Sf characterizes the Josephson 
current between the left and right wells. And 5"? is a displacement operator showing the 
imbalance between the wells. 
Let us introduce the notation 

Eo = ^{E, + E2) (12) 
and define the following combinations of the interaction matrix elements 

c.. - ^ (^r - v;r) ' ^.--2^r- (13) 

Also, we define the tunneling parameter 

n = E2-E^+Y^ dj . (14) 

Then Hamiltonian (6) reduces to the effective pseudospin representation 

H^EoN+^Y.AiJ - ^T.Sf + BijS^S] - E JijStS] , (15) 
in which N is the number of atoms in the lattice. 
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3 Evolution equations 



The Heisenberg equations of motion for operators (8) are 
Our aim is to study the evolution of the averages 

i i 

^^^j:<s!>, (17) 

i 

describing the temporal behavior of the tunneling intensity x = x{t), Josephson current 
y = y{t), and the imbalance variable z = z{t). 

Considering the statistical averaging in Eqs. (16), we need to choose a decouphng for 
the pair correlators of spin operators. The simplest option, which is usually accepted, would 
be the mean-field, or semiclassical, approximation, neglecting all pair correlations. Such 
a semiclassical approximation does not take into account quantum effects, due to particle 
collisions, which cause the attenuation of dynamical variables. 

To take into account such damping effects, we can employ the local-field approximation 
that is used in the theory of magnetic resonance [26,27] and in the calculations of the dielectric 
response functions [28]. The idea of this approximation is the existence of a kind of local 
equilibrium even in strongly nonequilibrium systems [29,30]. Particle collisions are considered 
as occurring in an effective local field of other particles [26]. Technically, for the system 
we consider here, the local-field approximation is realized as follows. First, one considers 
equilibrium mean-field solutions for Hamiltonian (15). Then, introducing the notation 

1 1 



we have 



N ^.-^ N ... 



n-2B<Sf> ^ ^(D 

— 2^^^^^n^ 



<5f>=0, >= ^^^^tanh(|;) , (19) 



where T is temperature and 



D = ^J{n-2B< Sf >)2 + 4J2 < S? >2 



The same form of solutions (19) is assumed to hold locally in time and space, when the 
system is in local equilibrium. This means that, in local equilibrium, the average variables 
(17) are of the form 

n-Bx , /A' 
Xt — — tanh I — 



A V2T 
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y* = , zt = ^ tanh ) , (20) 



A V2r 

in which 



The evolution equations for variables (17) are obtained by averaging Eqs. (16), decoupling 
the pair spin correlations, but including the damping, caused by particle collisions, so that 
the variables would attenuate to their local forms (20). This procedure yields the evolution 
equations 

— = Jyz- 72(x - xt) , 

dy 

— ^Qz-{J + B)xz - 72(2/ - yt) , 
dz 

— ^ -ny + Bxy - ji{z - zt) , (21) 

where the local fields are given in Eqs. (20). The attenuation parameters 71 and 72 can be 

calculated in the same way as for spin systems [27], being expressed through the particle 
interactions. In our case, these parameters arc of order J + B. 

Before solving Eqs. (21), it is convenient to reduce them to dimensionless notation. To 
this end, we define the dimensionless parameters 

uj = , h = (22) 

J + B' J + B ^ ' 

and notice that J/{J + B) — 1 — b. We also define the dimensionless function of time 

For the dimensionless damping parameter, we take 

7^^ = ^. (24) 
' J+B J+B ^ ' 

We consider the case of zero temperature and measure time in units oi 1/ {J + B). 

In this way, we come to the system of dimensionless equations describing the tunneling 
intensity, 

(1 - h)yz - 7(0; - Xt) , (25) 



Josephson current, 
and the well imbalance, 
with the local fields 



dt 

^ ^ {uj - x)z - -fy , (26) 

dz 

— ^{bx- uj)y - -f{z - zt) , (27) 

UJ — bx (1 — b)z , , 

= — , % = ^. (28) 
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where 

h = sjiuj - hxf + (1 - 6)2^2 _ (29) 

Equations (25) to (29) form the main system of evolution equations to be analyzed in what 
follows. We may notice that these equations are invariant under the change of signs of x, y, 
and LU. Fixing lu positive implies that the energy levels in a double well are enumerated so 
that Ei<E2. 



4 Nonlinear dynamics 

First, let us consider the existence and stability of attractors of Eqs. (25) to (27). We keep 
in mind that, by their definitions. 



a; > 



< 6 < 1 . 



(30) 



There are two fixed points that are either stable or unstable depending in the value of u. 
When the tunneling parameter is in the interval 



< a; < 1 , 



the stable fixed point is 



(31) 



(32) 



The characteristic exponents, defined as the eigenvalues of the Jacobian matrix, will be 
denoted by A. For the fixed point (32), we have 



Ai = -7 < 1) , 

and two other exponents are the solutions of the quadratic equation 
(1 - h)X^ + 7 (2 - 6 - cu') A + (l - cu') [(1 - hf + 7' 



. 



(33) 



(34) 



The general solutions of Eq. (34) are easily written, though they are rather cumbersome. 
For illustration, we represent the characteristic exponents for small 7 < 1. Then 



A 



2,3 



2(1 - h) 



7 ± i^eff , 



(35) 



with the effective frequency 



a;, 



eff 



(l-6)(l-a;2) 



1 



{b - cuYl^ 



8(1- 6)3(1 



(36) 



This shows that the fixed point (32), for small 7, is a stable focus. The value a; = 1 is a 
bifurcation point. When approaching this point, the characteristic exponents A2 and A3 are 



Ao~ 



2(1 -LO) 

(1 - bh 



+ 



;i-6)' + 7' 7^-4(1-6) 



(1 - 6)73 ' 
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X,^-^ + 2(1 - + [4(1 - bf + (3 + 6)f 1 ^^—^ , (37) 

for cj — s> 1 — 0. In the limit uj = 1. the fixed point (32) is neutral, since A2 = 0. When u 
grows larger than one, point (32) becomes unstable. 
For the tunneling parameter 

u;>l, (38) 

the stable fixed point is 

x*^l, y*^0, z*^0. (39) 
For the characteristic exponents, we have 

Ai = -7 {u;> 1) (40) 

and the equation 

A^+^-\-SA+'"-m<"-5+^^i=0 (41) 

UJ — b UJ — b 

defining A2 and A3. The latter for small 7 are 

2a; - 1 - 6 , . , 

with the effective frequency 

(1 - 6) V 



1 - 



8(cu-l)(cu-6)3 



(43) 



Recall that 6 < 1, hence uj > b. When uj approaches the bifurcation point from above, that 
is a; ^ 1 + 0, then 



A ~ (1-^)' + ^' iu, I) [(1-^)^ + 7^-27 . ..2 

(1-6)7 ^ ' (T^^W ^ ^ ^ ' 

A3^-7 + ^(c.-l) + ^^^^^(c.-l)^ (44) 

For UJ > 1, the fixed point (39) is stable. But at the bifurcation point uj = 1, the stationary 
solution (39) is neutral, since A2 = 0. And point (39) becomes unstable for a; < 1. 

An important observation is that the system of equations (25), (26), and (27) is struc- 
turally unstable if 7 = 0, since then Ai = and A2 and A3 are purely imaginary. Such 
a behavior is rather general for the cases not taking into account attenuation [31,32]. Al- 
lowance for the damping, due to particle collisions, makes the system structurally stable 
and substantially influences the overall dynamics. To emphasize the drastic change brought 
by the inclusion of the attenuation, we solve Eqs. (25) to (27) numerically and present the 
results in the figures. 

Figure 1 shows the system dynamics for the tunneling parameter in the region (31), 
when the stable fixed point is given by Eq. (32). The initial conditions correspond to atoms 
loaded into the lattice in a nonequilibrium state. If the attenuation is not taken into account. 
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the variables oscillate around their initial conditions. But giving due consideration for the 
attenuation results in the variables relaxing to the related stationary point. 

In Fig. 2, the parameters correspond to the stable fixed point (39). Again, at the initial 
time, atoms are loaded into the lattice in a nonequilibrium state. When there is no damping, 
all variables oscillate around their initial values. While, the presence of damping forces the 
variables to relax to their stationary solutions. 

Varying the parameter b in the admissible region < 6 < 1 does not influence much 
the system dynamics. Changing the tunneling parameter u inside the region (31) also does 
not yield considerable difference in the system behavior. When a; > 1, then its increase 
results in more pronounced oscillations of the dynamical variables. For comparison. Fig. 3 
presents the system dynamics, with the same initial conditions as in Fig. 2, but with a larger 
tunneling parameter. 

5 Conclusion 

An effective Hamiltonian for ultracold atoms in a double-well lattice is derived. Nonequi- 
librium situation is considered, when atoms are loaded into the lattice in a state far from 
a stationary one. As dynamical variables, we analyze the tunneling intensity, Josephson 
current, and the population displacement. The treatment is done in the local-field approxi- 
mation, taking into account the attenuation, caused by particle collisions. It is shown that 
the proper account of damping drastically changes the system dynamics, as compared to 
the case of no attenuation included. If the simple mean-field approximation, with no damp- 
ing, is used, the dynamic variables oscillate around their initial values. While taking into 
consideration the attenuation results in the relaxation of the dynamical variables to their 
stationary solutions. Depending on the value of the tunneling parameter, there exist two 
types of stationary solutions, whose main difference is either in the presence of a stationary 
atomic displacement in a double well, z]" 7^ 0, or in its absence, when the atoms in a double 
well are distributed symmetrically, Z2 — 0. 

An important message following from these results is that the so popular simple mean- 
field approximation should be used with caution. When considering nonlinear dynamics, such 
a simple approximation, with no account of attenuation, can lead to qualitatively incorrect 
conclusions. 

The results of the present paper are applicable to ultracold particles of different nature. 
These could be neutral atoms and molecules [7-11] or cold ions [33,34]. The basic requirement 
that the particles could be loaded into a double-well lattice. 

The double- well lattices possess a rich variety of properties, which makes them an attrac- 
tive candidate for different applications. A very important point is that the lattice properties 
can be regulated and controlled in a rather wide range. In addition to choosing the appropri- 
ate values of the lattice parameters, it is also feasible to realize their temporal modulation, 
similar to the case of usual optical lattices, whose parameters can be made time dependent 
by applying external fields [35-37], by moving [3,38], or just by shaking the lattice [39,40]. 
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Figure Captions 



Fig. 1. Dimensionless variables describing the tunneling intensity x{t), Josephson cur- 
rent y{t), and well population imbalance z{t) as functions of dimensionless time for a; = 0.1 
and b = 0.5. Initial conditions are Xq — 0.66, yo — 0.75, and zq = 0. The case of no attenua- 
tion (7 = 0) is shown by the dashed curve. The case with attenuation (7 = 1) is represented 
by the solid line. 

Fig. 2. Dimensionless variables x{t), y{t), and z{t) as functions of dimensionless time 
for uj = 1.5 and h — 0.5. Initial conditions are xq = 0.33, = 0.5, and zq — 0.8. The 
attenuation parameters are: 7 = (dashed line) and 7 = 1 (solid line). 

Fig. 3. Dimensionless tunneling intensity x{t), Josephson current y{t), and population 
imbalance (population displacement) z{t) as functions of dimensionless time for a; = 5 and 
b — 0.5. Initial conditions are xq = 0.33, yo = 0.5, and zq — 0.8. Behavior is shown for two 
cases: no attenuation, 7=0 (dashed line) and with damping, 7 = 1 (solid line). 
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Figure 1: Dimensionless variables describing the tunneling intensity x{t), Josephson current 
y{t), and well population imbalance z{t) as functions of dimensionless time for a; = 0.1 and 
b = 0.5. Initial conditions are xq = 0.66, yo — 0.75, and zq = 0. The case of no attenuation 
(7 = 0) is shown by the dashed curve. The case with attenuation (7 — 1) is represented by 
the solid line. 
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Figure 2: Dimensionless variables x{t), y{t), and z{t) as functions of dimensionless time 
for uj — 1.5 and b — 0.5. Initial conditions are xq — 0.33, yo = 0.5, and zq — 0.8. The 
attenuation parameters are: 7 = (dashed line) and 7 = 1 (solid line). 
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Figure 3: Dimensionless tunneling intensity x{t), Josephson current y{t), and population 
imbalance (population displacement) z{t) as functions of dimensionless time for a; = 5 and 
b — 0.5. Initial conditions arc xq = 0.33, yo = 0.5, and zq — 0.8. Behavior is shown for two 
cases: no attenuation, 7 = (dashed line) and with damping, 7 = 1 (solid line). 
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